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Abstract 

Terrestrial matter effects in neutrino propagation are T- invariant, so that any 
observed T violation when neutrinos pass through the Earth, such as an asym- 
metry between the transition probabilities P {i>n — > v e ) and P{y e — * i/„), would 
be a direct indication of T violation at the fundamental level. Matter effects 
do however modify the magnitudes of T- violating asymmetries, and it has long 
been known that resonant enhancement can lead to large effects for a range of 
plausible values of the relevant parameters. We note that the determinant of 
the commutator of the lepton mass matrices is invariant under matter effects 
and use this fact to derive a new expression for the T-violating asymmetries 
of neutrinos propagating through matter. We give some examples which could 
have physical relevance. 
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Evidence for oscillations of atmospheric neutrinos [jlj] indicates that at least muon 
and tau neutrinos participate in lepton mixing with large amplitudes. Other evidence 
of a deficit of detected solar neutrinos [Q] indicates that electron neutrinos must also 
participate in lepton mixing, and the most recent data || disfavour small mixing an- 
gles. The participation of all three species of neutrinos in lepton mixing raises the 
possibility of CP and T violations in neutrino oscillations, and, given the recent evi- 
dence of direct CP violation in the hadronic sector |4j], it would perhaps be surprising 
if such effects were not also manifested in the leptonic sector. The emergence of large 
mixing parameters in the leptonic sector offers the exciting prospect of potentially large 
CP- and T- violating asymmetries. One possibility |J which has not been excluded 
[[]] is that CP and T violations are maximal for neutrinos in vacuum. A number of 
other authors [[7| have explored the phenomenology of CP and T violation in neutrino 
oscillations for several different scenarios of lepton mass and mixing parameters. 

In both solar and atmospheric experiments, neutrinos can traverse a significant frac- 
tion of the Earth. As long-baseline accelerator and reactor experiments push to still 
longer baselines, the amount of material traversed by man-made neutrino beams also 
increases. This makes the understanding of matter effects in neutrino oscillations es- 
sential in order to determine the complete pattern of neutrino masses and mixings. The 
particle/anti-particle non-invariance of the matter term in the effective Hamiltonian 
for neutrino propagation in matter means that observation of a particle-antiparticle 
asymmetry does not necessarily indicate a fundamental CP- violation, although it is, 
in principle, possible to discriminate between the two sources of such effects using the 
observed distance- and energy-dependence ||. However, matter effects in neutrino 
propagation between two points at the surface of the Earth are T-symmetric, and 
any observed T violation, such as an asymmetry between the transition probabilities 
Piv^ — > u e ) and P[y e — > Vp) would certainly be a direct indication of T violation (and 
via CPT invariance, of CP violation) at the fundamental level. 

Despite the fact that matter effects in the Earth are T- invariant, they have a 
non-trivial effect on the signature of fundamental T violations of v and V separately, 
as well as on CP violations, through their influence on the effective neutrino mass 
and mixing parameters. In fact, it has long been known || that very large resonant 
enhancements of the T-violating asymmetries are possible in terrestrial matter. In 
this paper, we derive a simple new result for the effect of matter on the T- and 
CP- violating asymmetries, and explore some of the consequences in experimentally 
preferred scenarios of mass and mixing parameters. 

In the case of three generations of massive neutrinos, in a weak interaction basis 
which diagonalises the charged lepton mass matrix, Me = Dp, the neutrino mass ma- 
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trix, M v is in general, an arbitrary 3x3 matrix. The Hermitian square of the neutrino 
mass matrix, M u Ml, may be diagonalised to find its eigenvalues, and its eigenvec- 
tors form the columns of the lepton mixing matrix, U. It is well-known that under 
these circumstances, neutrinos propagating in vacuum undergo flavour oscillations, 
and furthermore, in general, these result in CP- and T-violating asymmetries. 

The CP- and T-violating asymmetries in the transition probabilities are given (for 
arbitrary mixing matrix) by the universal function 

P{y a -> up) - P(T7 a -+T7p) = P(v a -> vp) - P(up -> u a ) 

= 16Jsin(A 12 L/2)sin(A 23 L/2)sin(A3iL/2) (1) 

for any pair of flavour indices a, /3, where J is Jarlskog's mixing matrix-dependent 



invariant |10| and the Ajj = (Aj — Xj) are the three differences of eigenvalues of the 



Hamiltonian, H. In the vacuum case, it is sufficient to set 

H = M u Ml/2E (2) 

so that \i = mf/2E, where m; is the mass of the ith neutrino mass eigenstate (as 
usual, we number the eigenstates in increasing order of mass). We note that the 
parameter J has a maximal value of 1/(6a/3), and that the product of the three sine 
functions (the arguments are not independent) has a maximal value of 3\/3/8, which it 



takes when all three arguments are separated by 120° ||TTJ. The maximum magnitude 
of the product of three sines is controlled by the smallest of the three arguments, 
A12L/2, and we can consider 2/A12 as the reduced wavelength of the asymmetry. The 
asymmetry is observable only once this term has developed a significant phase, and 
if it is furthermore not averaged to zero by resolution effects. This fact limits the 
observability of CP and T violations in neutrino oscillations in vacuum to a window 
of parameter-space, as pointed-out in Ref. ||T2 . 



In the case that the neutrino beam passes through matter of uniform density, the 
Hamiltonian is modified to H' = H + AH, where 



AH = ±V2GN e 0. (3) 




The ± sign is to be taken as + for neutrino and — for anti-neutrino propagation, and 
makes explicit the particle/anti-particle asymmetry introduced by matter effects. The 
effect of matter is therefore to modify the mass eigenvalues and the mixing matrix 
elements, compared with their vacuum values. 

The matter electron density, N e , lies, for the Earth, in the range 0^iV e <6.2A^4cm _3 . 
The expression, Eq. (Q), for the T-violating asymmetries is still valid in the case of 
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propagation through matter of uniform density, except that the eigenvalues and the 
parameter J appearing in the expression are modified to their matter values A£ and J' 
respectively. For the CP-violating asymmetry between neutrinos and anti-neutrinos, 
the situation is made much more complicated by the particle/anti-particle asymmetry 
of the matter term in the Hamiltonian, and there are terms additional to that on the 
right-hand side of Eq. (|) [§. 

Jarlskog [ttDI has given an easy way to calculate the CP- and T- violation parameter, 
J, in terms of the mass matrices and their eigenvalues. This may be written: 

2A 12 A 23 A 3 i J = Im{Det[M, 2 , H]}/ (A e/Lt A^ r A Te ) (4) 

where the A^/ on the RHS refer to the differences between the squared charged lepton 
masses. This formula is valid in any weak basis, but it is useful for our purposes to 
evaluate the RHS in the weak basis which diagonalises the charged lepton mass matrix, 
in which it can be written: 

Im{Det[r^, #]}/(A eM A MT A Te ) = 2lm(H 12 H 23 H 31 ). (5) 

It is easy to see that use of the vacuum Hamiltonian, H, or the matter-modified one, 
H' in Eq. (|5p leaves the result invariant, as AH is diagonal in this basis. In fact, more 
generally, the commutator 

[DlH] = [DlH'} (6) 

is invariant to matter effects in this basis. Taking the determinant on both sides of 
Eq. (H) yields the physical result, valid in any weak basis, that Jarlskog's determinant 
of the commutator of the two lepton mass (-squared) matrices is invariant to matter 
effects, and, from Eq. (|) that: 

A 12 A 23 A 31 J = A' 12 A 

23^31 

J' (7) 

(where the primed quantities refer to the matter modified Hamiltonian, as above). The 
invariance, Eq. (0), is a principal result of this paper (seeming not to have appeared 
previously in the literature). One consequence is that it enables the phenomenology 
of T violation for neutrinos in matter to be stated in a new and transparent way. 

The T- violating asymmetries of Eq. (P can now be generalised to neutrino prop- 
agation through matter of uniform density in the following form: 

P\v a - up) - P'ivp - i/ a ) = leJ ^ 12 ^ 23 ^ 1 sin (A' 12 L/2) sin (A' 23 L/2) sin (A' 31 L/2) 

^12^23^31 

(8) 

where Eq. ([?]) has been used to rewrite J' in terms of its corresponding vacuum value, 
J, the vacuum masses, and the matter-modified masses, a form which is completely 
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independent of the matter-modified mixing matrix elements themselves. This result 
is exact, and valid for arbitrary vacuum Hamiltonian, ie. for arbitrary neutrino mass 
and mixing parameters. Eq. (^) shows clearly that in general, both the magnitude, 
and the wavelength of T-violating asymmetries are modified in matter compared with 
their vacuum values, and it makes explicit the correlation between these two effects. 
The modification of the magnitude is by a factor: 

J> _ A 12 A 23 A 31 

J A' 12 A' 23 A' 31 ' W 

which, in general, can be larger or smaller than unity. Provided that the asymmetry 
wavelength is not increased by matter effects beyond observability, potentially useful 
enhancements of the magnitude of T violation can occur in matter. 

In the form of Eq. (^), a number of results concerning T- violation for neutrinos 
propagating in matter can be seen rather easily. For example, it is manifest that even 
in matter, the T-violating asymmetries, as defined here, are independent of flavour. In 
the low-density limit, N e — > (Ay — > Ay for all the expression, Eq. (|j), reduces 
to the normal vacuum case, Eq. (ffl), as expected. Furthermore, in the small- L limit 
(ie. if the propagation distance is small compared with the scale of the shortest matter 
oscillation length, L « 2/A' 13 ), the equation reduces to the vacuum case, Eq. ([]]), 
as all the sine functions in Eq. (|§|) may be approximated by their arguments and all 
the primed variables cancel. For L — > therefore, there is no observable effect due 
to matter on T-asymmetries (nor indeed on any observable aspect of the oscillation 
phenomenon JT3J). In general, if J 7^ and there are no degeneracies in vacuum, 
Eq. (|7|) ensures that there are no degenerate eigenvalues in matter (A'^A'^A^ 7^ 0) 
and hence that TZ remains finite. 

We note that for an arbitrary Hermitian matrix, H, the combination of its eigen- 
values, Ai 2 A 23 A 3 i, is the square-root of the discriminant of its eigenvalue equation: 

A 12 A 23 A 31 = V-4S 3 + S 2 T 2 - 27V 2 + 18VTS - 4X?T 3 (10) 

where the invariants T, S and D are respectively the trace, the sum of the principle 
minors and the determinant of H. This can be further simplified by noting that we 
have the freedom to add an arbitrary multiple of the unit matrix to H without altering 
its discriminant, so that we can choose either T = or one of the eigenvalues, A; = 0. 
So, for any given model of the neutrino masses and mixing angles in vacuum, the 
magnitudes of the T-violating asymmetries in matter can, in principle, be calculated 
directly from the vacuum parameters of the model and the (simplified) invariants of the 
matter-modified Hamiltonian, H', thereby obviating the need to diagonalise explicitly 
H' to find the matter-modified mixing angles. 
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In the general case, application of Eq. flUf ) to the matter-modified Hamiltonian to 
calculate the denominator of Eq. (^) yields the square- root of a quartic function of N e , 
with coefficients which are functions of the vacuum mixing parameters (masses and 
mixing angles) and the neutrino energy. The presence of this quartic function means 
that TZ has either one or two resonant maxima, as a function of the matter density. 
These correspond to the well-known matter resonances, occuring approximately where 
pairs of matter- modified neutrino masses are most nearly degenerate [O]. Values of 
TZ can be arbitrarily large, if J is sufficienctly small (clearly, J' = RJ cannot exceed 
1/6^3). For very large values of N e , the quartic term dominates, and T asymmetries 
are suppressed by ~ 1 /N%. 

In the remainder of this paper, we explore some specific examples which could have 
physical relevance. In order to ensure that our discussion is relevant to experiment, 
we will restrict our considerations to regions of parameter space which are not ex- 
cluded by present, corroborated neutrino experiments, such as the Super Kamiokande 
atmospheric neutrino data |1[] and the Super Kamiokande and Gallium solar neutrino 
data 0. We will use for the lepton mixing matrix, U, a conventional form in which 
there are three real mixing angles and one complex phase such that: U e s = sin# 13 , 
U e 2 = sin 9 12 cos 6*i3, = sin ^23 cos ^i3e J<5 and all the other elements are fixed by 
unitarity. 

A priori, our preferred scenario for neutrino oscillations is the threefold maximal, 
or tri- maximal scheme ||, which still provides a broadly consistent account of all 
the corroborated sets of experimental data on neutrino oscillations. In this scheme, 
all the elements of the lepton mixing matrix have equal moduli of magnitude 1 / v3, 
the vacuum value of J takes its maximal value, l/(6v^3), and all observed neutrino 
disappearance data are the result of neutrino oscillations governed by the scale of the 
larger neutrino mass-squared difference, <5m 2 3 ~ 1CT 3 eV 2 , with 5m\ 2 unresolved even 
by the solar data. This scheme is summarised by 12 = #23 = 45°, #13 = sin" 1 (l/y/3) 
and 5 = 90°. In this case, TZ is always less than unity in matter, so that asymmetries 
are always suppressed. 

A viable, and even experimentally favoured, alternative to tri-maximal mixing 
effectively factorises the atmospheric and solar scales by setting f/ e3 ~ in the mixing 
matrix fl4| . The atmospheric data then require #23 — 45°, and an energy-independent 
solar supression of 1/2 is obtained by setting #12 = 45° in the original bi-maximal 



scheme [pf] . We have ourselves proposed (see also |15|]) a particular variant of the 
bi-maximal scheme with ?7 e2 = — U T 2 = 1/3 which simulates very effectively tri- 
maximal mixing with a solar supression of 5/9. In bi-maximal-type schemes, Jcp = 
and there are no fundamental CP- or T-violating asymmetries involving neutrinos, in 
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vacuum or in matter. There are still fake vjv asymmetries due to matter effects, but 
these are of little fundamental interest. 

For illustrative purposes, we have explored in this paper an ansatz which interpo- 
lates smoothly between tri-maximal and bi-maximal mixing: 



with S = 90°. We have assumed that <5m 2 3 ~ 10 3 eV 2 and have allowed 5m\ 2 to 
remain variable. This scheme has the following properties: 

• and v T are treated democratically 

• v\ and U2 are treated democratically 

• in the limit sin 2 13 — > 1/3, we obtain tri-maximal mixing 

• in the limit sin 2 #13 — ► 0, we obtain bi-maximal mixing 

• Jcp — (1/4) sin#i 3 cos 2 0i 3 varies between its minimal and maximal values, as 
we move from bi-maximal to tri-maximal mixing. 

For arbitrary sin# 13 (0 < sin# 13 < ^=), there is always a resonance where T- violating 
asymmetries are maximised and this can be probed by choosing the neutrino energy, 
E, or the density of matter traversed, appropriately, as long as there is sufficient 
pathlength for the asymmetry to develop. The values of neutrino energy and/or matter 
density at resonance, and the maximum magnitude of the asymmetry there depend 
on sin #13 and the vacuum masses. 

Figs, la and lb show some examples of the maximum magnitude of the T-violating 
asymmetry defined in Eq. (H), as a function of the matter density in units of Na cm -3 
for several values of sin# 13 in this ansatz. The two figures differ in terms of the 
hierarchy of vacuum mass values used, and illustrate cases with one and with two 
resonant densities respectively. All the curves are for a fixed neutrino energy, E — 1.5 
GeV, which turns out to be the most suitable energy scale to maximise these effects in 
terrestrial matter within this model (see Eq. ([IB])). The Earth's mantle has an electron 
density of approximately 2Na cm -3 and the core, roughly QNa cm -3 . The negative 
density parts of the curves correspond to anti-neutrinos propagating in normal matter 
with electron density \N e \. We note that enhancements for neutrinos, for example, 
are typically compensated by suppressions for anti- neutrinos, and/or for neutrinos 
at different energies. As shown in Ref. ||, the electron neutrino begins to decouple 



U 



I ^= cos 6»i3 -^cos6» 13 

-i(l + sin# 13 e i5 ) i(l-sin# 13 e i5 ) 
v |(1 -sin# 13 e i5 ) -f(l + sin# 13 e i5 ) 



^cos0 13 e i<5 
^ cos 0i3 e iS j 



sin 6>i3 \ 



(11) 
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completely in the high energy limit, and T asymmetries are suppressed accordingly. 
The magnitudes of the asymmetries in vacuum can be read from the curves at the zero 
density point. Several other features are typical of the generic MSW-like matter effect, 
eg., the resonance gets sharper and the asymmetry maximum gets closer to unity as 
the vacuum mixing angle decreases. 

It is interesting to consider under what circumstances the asymmetry becomes 
maximal in this ansatz. Even in the limit sin #13 — > (sin #13 7^ 0), where T- and 
CP-violating asymmetries in vacuum become arbitrarily small, they can reach their 
maximum possible value of unity in matter. We find that this can be achieved if 

x 2 a \Z2-cos0i3 2 

bm 12 = 2 cos 6>i3— — — bm 13 (12) 
1 + sir 6>i3 

and then the matter density and neutrino energy at resonance are related by: 

y/2 l + sin 2 fl 13 13/ V ; 

Under the above two conditions, threefold maximal mixing is achieved at resonance. 
Eq. fll2"|) implies that for such mixing at resonance, the two Am 2 values must be of 
the same order of magnitude, which is probably not ruled out by the present data. 

Fig. 2 shows an example of how one of the asymmetries in Fig. 1 develops with 
propagation distance, in units of 1000 km. It is plotted for neutrinos in matter with 
a constant density of 1.9Na cm -3 , ie. close to the density of the Earth's mantle. The 
matter-enhanced asymmetry is compared with the vacuum asymmetry and the corre- 
sponding matter-enhanced asymmetry for anti-neutrinos. The enhanced asymmetry 
does not develop to the 100% level within the distance scale of the Earth's diameter in 
this case, (it would, given a longer pathlength), although there are cases where it does. 
It does however exceed considerably the vacuum asymmetry, along almost the whole 
trajectory within the Earth, and may be observable when the vacuum asymmetry 
would not be. 

While Eq. ( |T2] ) may not be satisfied exactly in nature, there is a significant range of 
parameter space over which large enhancements are possible. We have found the con- 
ditions under which the enhancements are maximal, and here, at least for some values 
of the parameters, matter can induce a situation where the mixing matrix is arbitrarily 
close to the tri-maximal form. Such a scenario might well be obtained in nature, and 
we recommend that large T- and CP-violating asymmetries be searched for in future 
long-baseline neutrino experiments. It may even be possible to site experiments to 
exploit the matter effects in the Earth so as to maximise the asymmetries. 
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Figure Captions 



Figure 1. Examples of resonant enhancement and suppression of T violation for neu- 
trino oscillations in matter using the ansatz of Eq. ([0]). The maximum magnitude of 
the T- violating asymmetry is plotted as a function of matter density and sin# 13 , for 
neutrino energy, E = 1.5 GeV, with a). <5m 2 3 = 1.0 x 1CT 3 eV 2 and 5m\ 2 = 0.7 x 10~ 3 
eV 2 ; b). <5m 2 3 = 1.0 x 10~ 3 eV 2 and 5m\ 2 = 0.2 x 10~ 3 eV 2 . In each case, the point 
marked "o" represents the vacuum value of the asymmetry for the same vacuum mass 
parameters. The Earth's mantle has an electron density of approximately 2Na cm -3 
and the core, roughly 6Na cm -3 . NB. The density scale can be converted to a scale 
of energy in GeV for neutrinos propagating in the Earth's mantle by multiplying the 
numbers by ~ 0.75. 

Figure 2. An example of resonant enhancement of T violation for neutrino oscillations 
in matter. The solid curve shows the T-violating asymmetry as a function of propa- 
gation distance, for neutrinos of energy, E = 1.5 GeV, with 5m 2 3 = 1.0 x 10~ 3 eV 2 , 
5m\ 2 = 0.7 x 10~ 3 eV 2 and sin #13 = 0.1. The matter electron density in this example 
is 1. 97X^4 cm -3 . The dashed line shows the same quantity in vacuum for the same 
vacuum parameters, and the dotted line shows the same quantity for anti-neutrinos. 



11 



Max. Amplitude of Asymmetry 
p p p p 
o ho *-f^ en bo — >■ ho 
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